007 Nature Publishing Group _http://www.nature.com/naturegenetics 


CORRESPONDENCE 


except the group of university students (P = 
0.028), were in Hardy-Weinberg equilibrium 
(by a Markov chain method in GENEPOP). 
Furthermore, we compared the genotypes 
among individuals with SARS with differ- 
ent prognoses. If L-SIGN homozygosity is 
a protective factor against infection, it may 
also be associated with better prognosis after 
acquiring the infection. Therefore, we also 
examined whether homozygotes had a bet- 
ter prognosis by classifying individuals with 
SARS who had an uneventful recovery versus 
those who had severe disease and wereadmit- 
ted to the intensive care unit for mechanical 
ventilation support (an approach similar 
to that reported previously®). However, we 
did not detect any significant association 
(P =0.9, Supplementary Table 1 online). 

Sample size is the main limitation of both 
studies. However, these two samples already 
represent the few ‘large’ collections of indi- 
viduals with SARS available for genetic study. 
To estimate the size of an overall effect, we 
performed a meta-analysis of thetwo datasets 
together by the Mantel-Haenszel test using 
control groupsin Hardy-Weinberg equilibrium 
(two groups of controls in this study and ran- 
dom controlsfrom Chan et al.; total n = 1,468; 
462 affected individuals and 1,006 controls). 
The combined odds ratio was not significant 
(combined OR = 0.84; 95% confidence inter- 
val: 0.66- 1.06, P = 0.14). 

The difference in the results between the 
two studies was basically accounted for by 


a difference in the homozygote proportions 
in the controls (45.4% in this study versus 
55.0% in Chan et al.), while the homozy- 
gote proportions among individuals with 
SARS are almost identical (46.9% here ver- 
sus 46.3% in Chan et al.). The reason for the 
discrepancy in the homozygote proportions 
in the ‘control’ groups is not clear. However, 
a subpopulation difference in allelic and 
genotypic frequencies exists between north- 
ern and southern Chinese. The seven-repeat 
allele was more prevalent in the Beijing 
sample (0.7 in Beijing versus 0.64 in Hong 
Kong; P =0.05), which also largely accounted 
for the higher proportion of homozygotes 
(55.7% in Beijing versus 46.0% in Hong 
Kong; P = 0.02). Unrecognized subpopula- 
tion structure may confound genetic asso- 
ciation studies. Results in the study by Chan 
et al. suggested that this confounding factor 
might be present. There were three groups 
of controls, including two groups of hospital 
controls (health care workers who worked in 
SARS wards and affected individuals attend- 
ing various outpatient clinics) and a group of 
blood donor controls. Interestingly, genotype 
distributions from both groups of hospital- 
based controls deviated significantly (P < 
0.0001) or marginally (P =0.05) from H ardy- 
Weinberg equilibrium. 

In addition, other yet-unknown mecha- 
nisms (such as alternative splicing of theneck 
region, which could interfere with formation 
of homotetramers among homozygotes) may 


account for the discrepancy between the two 
studies. Replication isan important approach 
to verify any significant genetic association 
findings®’, and additional association studies 
are required to establish the putative protec- 
tive effect of L-SIGN homozygosity against 
SARS or other infections. 
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© Lack of support for an association between CLEC4M 


coronavirus infection 


To the Editor: 

Inthe] anuary 2006 issue, Chan et al. reoorted 
a significant association between severe acute 
respiratory syndrome (SARS) and a variable 
number of tandem repeats (VNTR) polymor- 
phism in exon 4 of CLEC4M in a collection 
of individuals from Hong Kong. CLEC4M 
encodes L-SIGN (‘liver/lymph node specific 
ICAM -3 grabbing nonintegrin’), which serves 
as areceptor for many viruses, including SARS 
coronavirus (CoV)2. Individuals homozygous 
for CLEC4M tandem repeats were reported to 
beless susceptibleto SARS CoV infection. The 
authorsalso showed that cells homozygous for 
CLEC4M repeats had a higher binding capac- 
ity for SARS CoV, higher proteasome-depen- 
dent viral degradation and alower capacity for 


trans infection. Thus, both genetic and func- 
tional studies suggested that homozygosity 
for CLEC4M was associated with protection 
against SARS CoV infection. 
Itisimportant to bear in mind that associa- 
tion studies requirereplication in independent 
populations? Wetherefore attempted to repli- 
cate the findings of Chan et al. by genotyping 
the VNTR polymorphism in three additional 
collections of case-control samples from 
northern China: (i) the ‘Beijing community 
population’, consisting of 339 individuals with 
SARS and 227 random controls recruited from 
the community‘; (ii) the ‘Beijing health care 
worker (HCW) population’ consisting of 42 
health careworkers infected with SARS during 
the course of hospital duty and 40 health care 


@homozygosity and protection against SARS 


workers who had worked in SARS wards but 
remained free of disease and wereconfirmed to 
beseronegative for SARS? and (iii) the’Tianjin 
population’, consisting of 60 individuals with 
SARS and 129 disease-free controls (including 
85 random controls and 44 health care work- 
ers)®, Thethree collections of case-control sam- 
ples and their ascertainment criteria have been 
described in detail previously (Supplementary 
Methods online)*®. All groups except the 
individuals with SARS from the Beijing com- 
munity werein Hardy-Weinberg equilibrium. 
We found no significant differences in allde, 
genotype and homozygote or heterozygote 
frequencies between affected individuals and 
controls in the three populations (Table 1 and 
Supplementary Table 1 online). Early reports 
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Table 1 Homozygote or heterozygote frequencies of the CLEC4M polymorphism in sam- 


ples from northem China 


VNTR polymorphism? 


Heterozygotes 


Homozygotes 


Beijing community population 
All cases (n = 339) 

Random controls (n = 227) 
OR (95% c.i.) 

P value 


Cases with comorbid conditions (n =52) 
Random controls (n = 227) 

OR (95% c.i.) 

P value 


Cases without comorbid conditions (n = 287) 
Random controls (n = 227) 

OR (95% c.i.) 

P value 


Severe cases (n = 19) 
Mild cases (n = 268) 
OR (95% c.i.) 

P value 


Beijing HCW population 
HCW cases (n =42) 
HCW controls (n = 40) 
OR (95% c.i.) 

P value 


Tianjin population 
All cases (n = 60) 

All controls (n =129) 
OR (95% c.i.) 

P value 


200 (59.0%) 
117 (51.5%) 


139 (41.0%) 
110 (48.5%) 
1.29 (0.89-1.87) 
0.19 


20 (38.5%) 32 (61.5%) 
110 (48.5%) 117 (51.5%) 
1.36 (0.60-3.07) 

0.46 


119 (41.5%) 
110 (48.5%) 
1.26 (0.86-1.85) 
0.23 


168 (58.5%) 
117 (51.5%) 


8 (42.1%) 
111 (41.4%) 
1. 00 (0.38-2.59) 


11 (57.9%) 
157 (58.6%) 


1.00 
14 (33.3%) 28 (66.7%) 
18 (45.0%) 22 (55.0%) 
1.51 (0.58-3.99) 
0.40 
33 (55.0%) 27 (45.0%) 


72 (55.8%) 57 (44.2%) 
1.18 (0.60-2.34) 
0.63 


@ 


OR, odds ratio; c.i., confidence interval; VNTR, variable number tandem repeat; HCW, health care worker. 
*The heterozygotes are used as the reference group, and all ORs and P values are adjusted for age and gender. 


Primer sequences used for genotyping are listed in Supplementary Table 2 online. 


have shown that somenongenetic factors, such 
as comorbid conditions (including diabetes 
mellitus, hypertension, heart disease, tuber- 
culosis, asthma and malignancy), are risk fac- 
tors for the development of SARS”8 and may 
confound the contribution of genetic factors 
to thisdisorder. H owever, after stratification by 
comorbid conditions, theassociation remained 
nonsignificant in our Beijing community 
population (Table 1). Of the 287 affected indi- 
viduals without comorbid conditions, 19 were 
individuals with severe SARS who wereidenti- 
fied by their admission to intensive care units 
or by their death, and the remaining 268 were 
individuals with mildly symptomatic SARS. 
To account for this, we assessed whether there 
was an association between homozygosity for 


CLEC4M tandem repeats and severity of SARS, 
but we found none (Table 1). This result may 
be due to the limited number of individuals 
with severeSARS in thecurrent study and will 
require confirmation in additional studies. 
Collectively, theresultsin our threecollections 
of case-control samples from northern China 
are not supportive of thefindings of significant 
association between the VNTR polymorphism 
and SARS risk reported by Chan et al.1. There 
are several possible reasons for the inconsis- 
tent results. First, inadequate power may bean 
explanation of our negativeresults. However, the 
sample size in our Beijing community popula- 
tion (with 339 affected individualsand 227 con- 
trols) is approximately similar to that used by 
Chan et al. (with 285 affected individuals and 
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380 random controls), and this sample sizehad 
power >0.94 to replicate the effects by Chan et 
al. (calculated by the genetic power calculator?). 
Additionally, just by using our Beijing commu- 
nity sample set, we have successfully confirmed 
the positiveassociations between the mannose- 
binding lectin polymorphisms and SARS risk* 
that were observed previously in a case-control 
sample from Hong Kong2, Furthermore, the 
consistency of the negative associations in our 
Beijing HCW population and our Tianjin popu- 
lation strengthens our results. 

Second, there may be a small, population- 
specific difference in the contribution of 
CLEC4M polymorphism to SARS susceptibil- 
ity. This might occur if there were population 
differences in linkage disequilibrium pattern 
or allele frequencies of CLEC4M . Indeed, the 
homozygote or heterozygote frequencies in 
individuals with SARS in our Beijing com- 
munity or Beijing HCW populations were 
significantly different from those reported for 
affected individuals in Hong Kongcommunity 
or HongKongHCW populations (P = 0.0066 
and P = 0.017, respectively; Supplementary 
Table 1). Additionally, the difference in allele 
frequency between Hong Kong outpatient 
controls and Tianjin random controls and 
the differences in homozygote or heterozy- 
gote frequencies between Hong Kong HCW 
controls and Tianjin HCW controls were also 
significant (P = 0.0064 and P = 0.017, respec- 
tively; Supplementary Table 1). Furthermore, 
there was also a significant difference in the 
genotype and homozygote or heterozygote 
frequency between the Chinese population 
and Europeans!. 

Another possibility isthat Chan et a/. unwit- 
tingly neglected to account for the potential 
confounding factors that may distort the con- 
tribution of CLEC4M VNTR polymorphism to 
SARS susceptibility. Although Chan et al. took 
stringent precautions to stratify their samples 
by ‘health care worker’, population samples 
may still differ by many other factors that 
depend on setting and context of recruitment, 
such as age of the subjects at SARS onset, sex 
of the subjects and any comorbid conditions. 
Unfortunately, Chan et al. do not provide data 
with regard to such important information on 
confounding factors. Last, the initial findings of 
Chan et al, may not represent real associations 
and might be false positives. In genetic associa- 
tion studiesof common diseases, thereisa very 
low prior probability of detecting a true asso- 
ciation result when accounting for statistical 
adjustment for multiple comparisons. Indeed, 
the genetic associations presented by Chan 
et al. were marginally statistically significant 
(P = 0.027, 0.045 and 0.031, when compar- 
ing all SARS samples to random controls, 
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community SARS to outpatient controls and 
HCW SARS to HCW controls, respectively). 

Thus, we did not find any significant differ- 
ences in allele, genotype and homozygote or 
heterozygote frequencies between cases and 
controlsin our threeindependent populations 
of northern Chinese. Although the biological 
plausibility of L-SIGN and the functional evi- 
denceof theVNTR polymorphism in the origi- 
nal report remain interesting, we urge that the 
association between CLEC4M polymorphism 
and SARS be investigated in other subpopu- 
lations of ethnic Chinese origin (for example, 
Taiwanese or Guangdong Chinese) or in those 
of different ancestry, such as Europeans. 
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Chan etal. reply: 

Tang et al. and Zhi et al. report that in 
contrast to our findings?, they were unable 
to find association between homozygosity 
or heterozygosity of theCLEC 4M (L-SIGN) 
exon 4 tandem repeat polymorphism and 
SARS CoV susceptibility in the Chinese 
population. Their data cannot conclusively 
negate our findings for the reasons below. 

Weagree with Tang et al. that the 
difference between their analysis and ours 
is largely accounted for by the difference 
in the percentage of homozygosity in 
controls. Tang's controls included neonatal 
cord blood, healthy elderly individuals 


Je) aged >70 years, local university students 


and control samples collected in Beijing. 
Thehomozygosity and heterozygosity 
distribution of their Beijing controlsis very 
similar to that of our Hong Kong random 
controls. H owever, the former three control 
groups collected in Hong Kong are poorly 
matched by age for comparison with 
individuals with SARS and our random 
controls (Supplementary Table 1 online), 
and such an age distribution clearly excludes 
the population aged 25-69 years, which 
makes up the largest proportion (>70%) 

of SARS-infected individuals during the 
Hong Kong outbreak in 2003 (ref. 2). Age 
and gender are well-known confounding 
factors in any case-control study, and for 
SARS, individuals aged >65 years and <18 
years are actually associated with a lower 
risk3, Thus, itis surprising that Tang et 

al. chose these specific control groups at 
extreme age ranges, given that there areno 
obvious obstacles for recruiting appropriate 


age-matched controls in Hong Kong. They 
also failed to perform logistic regression 
accounting for age differences, which should 
have been incorporated in their statistical 
analysis. T he age and gender distribution 

of the SARS and control populationsin our 
study, in contrast, is more properly matched, 
and by logistic regression, our results 
remain statistically significant (Table 1 and 
Supplementary Table 1). Our results remain 
consistent after accounting for comorbid 
conditions (Table 1). 

Tang et al. also used control groups 
differing widely in age to determine if age 
had any effect on genotype frequencies. 
However, their data per se already showed 
asignificant difference in overall genotype 
distribution, when comparing neonates 
versus elderly individuals versus university 
students (P = 0.029, y? test). The difference 
in genotype between neonates and university 
students is also significant (P = 0.009 
by CLUMP), suggesting that age-related 
selection may exist for CLEC4M genotypes. 
Indeed, age-dependent variation of allele 
and genotype frequencies has been reported 
for other genes*5, 

Tanget al. also claim to show a 
subpopulation differencein allele and 
genotype frequencies between northern 
and southern Chinese. T hey note that “the 
seven-repeat allele was more prevalent in 
the Beijing sample (0.7 in Beijing versus 
0.64 in HongKong, P = 0.05), which also 
largely accounted for the higher proportion 
of homozygotes (55.7% in Beijing versus 
46.0% in Hong Kong, P = 0.02).” It seems 
these quoted H ong Kong figures refer to 


frequencies from their cord blood group 
alone. Such a comparison is difficult to 
justify, again because age is seriously 
mismatched. It should be noted that thereis 
no differencein either the seven-repeat allele 
frequency or the proportion of homozygotes 
between the H ong Kong random controls 
of Chan etal.) and the Beijing controls of 
Tang et al. (Supplementary Table 2 online). 
Barreiro et al. investigated the CLEC4M 
homozygote and heterozygote distribution 
of different ethnic groups and reported 

that the proportion of homozygotes in 

East Asians, consisting predominantly of 
Chinese, was 53% (Supplementary Table 

3 online), a figure similar to the Hong 
Kongrandom controls of Chan et al. and 
the Beijing controls of Tanget al. Thereis 
also no significant differencein allele or 
genotype frequencies of these two control 
populations from the East Asian samples 

of Barreiro etal.’ (Supplementary Table 
2).When these results are taken together, 
little subpopulation structure is observed 
for CLEC4M homozygote and heterozygote 
distribution between the northern and 
southern Chinese populations. 

Finally, given the apparent absence of 
subpopulation structure for CLEC4M 
homozygote and heterozygote distribution 
in the Chinese population, we performed 
a meta-analysis of our data set and that 
of Tang et al. by the Mantel-H aenszel test 
usingall control groups that arein Hardy- 
Weinberg equilibrium: our random controls 
and outpatient controls (n = 670), and 
cord blood, healthy elderly individuals 
and Beijing controls (n = 827) of Tang 
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